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Abstract 

The conventional method to determine the cosmic curvature is to measure 
the total mass density fitot- Unfortunately the observational fitot is closely near 
the critical value 1. The computation of this paper shows that f2tot ~ 1 is an 
inevitable result for the young universe independent of the spatial topology. So 
the mass density is not a good criterion to determine the cosmic curvature. In this 
paper, we derive a new criterion based on the galactic distribution with respect 
to redshift z, which only depends on the cosmological principle and geometry. 
The different type of spatial topology will give different results, then the case 
can be definitely determined. 



PACS numbers: 98.80.-k, 98.65.-r, 98.80.Es, 98.80.Jk 



Key Words: cosmic curvature, large scale structure, galactic 
distribution, galaxy count 

* email: yqgu@fudan.edu.cn 

t email: Maxim. Khlopov@romal .infn.it 



1 



1 Introduction 



Whether the space of our universe is open, flat or closed is one of the most fundamental 
problems in cosmology. The solutions to some other important problems such as the 
property of dark matter, the evolution of the universe, are closely related to the spatial 
topology of the universe. Almost all used methods to determine the cosmic curvature 
are based on computing the total mass density f2 to t- However, the most authoritative 
empirical data show that the total mass density f2 to t = 1.02±0.02[l]-[4], which is closely 
near the critical density 1, so we can not get a definite solution to this fundamental 
problem. 

Recent years, lots of attention have been attracted to the resolution of this problem [5]- 
[15]. In [7], by analyzing the large-scale correlation function of red galaxies from 
the Sloan Digital Sky Survey(SDSS) and combining with the cosmic microwave back- 
ground(CMB) acoustic scale, we arrive at Qk — —0.010 ± 0.009 under the assumption 
of w = —1. In [8], by measured the large-scale power spectrum P(k) for the luminous 
red galaxies from SDSS, the authors improve the evidence for spatial flatness, and 
sharpen the curvature constraint to Q tot = 1.003 ± 0.010. 

In [9, 10], using the full dataset of CMB and large scale structure(LSS), and by 
Markov Chain Monte Carlo global fit, the authors find the best fit value of \Qk\ < 0.015 
for some dark energy models and \Qk\ > 0.06 are excluded. The ACDM model is 
consistent with all the data, and a flat universe is preferred. 

From the calculation of the above papers, we find that the results are strongly and 
nonlinearly correlated with the equation of state w(z). In [14], the authors solved w(z) f 
and w(z)d, two equations of state for the same model, respectively from the Friedmann 
equation and the luminosity distance equation. The numerical results show that the 
two functions w(z) unstably depend on Qk and depend on z towards opposite trends, 
which lead to great errors for adequately large z. Similar conclusion also confirmed 
in [15], where the authors find bounds on cosmic curvature are less stringent if dark 
energy density is a free function of cosmic time. 

However, in our point of view, Q t ot or VLk may be not the best criterion to determine 
the cosmic curvature. At first, Q to t ~ 1 is an inevitable results for the young universe 
as shown bellow. Secondly, since the Friedmann equation is a dynamic equation of a(t), 
it is inadequate to analyze this equation as algebraic equation by introducing Hubble's 
parameter H . Thus looking for a more effective method and criterion to determine the 
spatial type is necessary. 
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Hubble once realized to measure the spatial curvature by counting galaxies[16], but 
the corresponding equation derived is very complicated, which not only depends on z 
but also on H and Q t ot- This situation is caused by the calculation, which is too direct 
to simplify the model, so that the idea was not developed for practical application. 

However, the geometrical problem may be more naturally and effectively solved by 
geometrical method. In this paper we develop the Hubble's idea by designing a simple 
reference function of the galactic distribution in flat space, and then compare it with 
the practical counting data. The result of comparison may be able to give a definite 
answer to the fundamental question. 

2 Analysis and method 

2.1 Flatness problem of the space 

In mean sense, the universe is highly homogeneous and isotropic, and the metric is 
described by Friedmann-Robertson-Walker(FRW) metric. The corresponding line ele- 
ment is given by 



where K — 1, and —1 correspond to the closed, flat and open universe respectively. 
Some literatures intended to replace the constant K by some continuous function K(r) 
or K{t). In fact, such replacement violate the symmetry of the spacetime. The cosmic 
radius of curvature is determined by a{t) rather than by K. The spatial topology of 
the homogeneous and isotropic spacetime has only the three cases[17]-[20], and the 
metric must be equivalent to (2.1). In this paper, we adopt the conformal coordinate 
system, because in this form a(t) can be expressed as Taylor series near a — > 0. The 
corresponding metric becomes 




(2.1) 




(2.2) 




(2.3) 



1 sinhr if K — — 1. 
The critical density p c and the Hubble's parameter H are defined by 




(2.4) 
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By Friedmann equation[16], we have 



n K = K 




vt m + i^a _ i — "tot _ i) 



(2.5) 



where Vt m = &>a = ^pp- Theoretically, by (2.5) we have the following judgment: 
Q tot > 1, Q tot = 1 and Q tot < 1 correspond to the closed, flat and open universe 
respectively. The following analysis shows that, this judgement is not sharp for a 
young universe. 

The original equation of (2.5) is the dynamic equation 



For most matter models [21], such as nonlinear spinors[22, 23, 24] and Casimir effect, 
when a — > 0, we have the main part of density as 



where p is a constant, and \a\ <C a is a slowly varying function of a, which can be 
treated as a constant for the following calculation. For homogeneous scalar field models 
such as quintessence and phantom field [25], we can also get the corresponding p(a) 
via transforming ^ = ^a' and then solving the differential equation. However, we 
need not to solve the specific model for the present purpose, because we only need the 
assumption that the cosmic equation of state p(a) exists and has singularity not worse 
than (2.7). There are a lot of works on constructing complicated equation of state w(z) 
and fitting the empirical data. It seems unnecessary, because we can hardly find out 
the effective w(z) [14]. 

We take (2.7) with constant a as model to show the properties of a(t) and Qk 
for small t. We take R = ^y^Po = 1 as length unit, which is the mean scale factor. 
Substituting (2.7) into (2.6) we get 



a 



.12 




(2.6) 




(2.7) 



a 



2(a -a)- Ka 2 + Aa 4 , 



(2.8) 



where A = |Ai? 2 is dimensionless constant. Choosing suitable starting time t = 0, for 
small t the solution of (2.8) can be expressed by 




( K ( K 2 \ 3A 2 \K 3 \ 4 ^. 6N 
(-24 + ( Xa + 24 ) a ° ~ T a ° + ^ al ) 1 + 
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Figure 1: The Maximum Deviation of Mass Desity Caused by Space Type 

where < |ao| ~ \<j\ <C 1 is the root of 2 (a — a) — Ka 2 + Aa 4 = near zero. For the Big 
Bang model we have ao = 0. (2.9) shows how the parameters K, a, A and a influence 
the scale factor. 

As computed in [22, 23], we have the following conclusions. (I) Since A has effects 
only on the cosmic scale, we have |A| < 1. (II) The present comving time t a ~ 18° = 
0.314(rad). In this paper, we are only concerned with the behavior of the universe after 
the galaxies formed, so we have the estimation for maximum of as follows 

i^IHH) 2 ^' - - 1 - K :: + 2 H - ° (2 - io) 

The typical values are displayed in Fig.(l), which visually tells us VLk is not a good 
criterion to determine the cosmic curvature. If ao < 0, then flx{t) = has solution 
^singular — ^> which inevitably leads to the so called "fine-tuning mass density", thus 
the dynamic process can not be understood with static equation. 

2.2 Galactic distribution and reference function 

Now we design the new method to determine the cosmic curvature K, which only 
depends on the redshift distribution of galaxies at a few selected point zj.. The cosmo- 
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logical redshift 



z = 



A — A a(t) 



1, 



(2.11) 



Ao a(t ) 

where t = t + r is the present time, and to is the time when the photon was emitted, r 
is the radial coordinate of the light source. Denote the number of all galaxies in a ball 
with radius coordinate r by N(r), we consider N(r) in the domain [r — |Ar, r + |Ar], 
or correspondingly in [z — \5z, z + \5z]. The increment of redshift is given by 



a(t ) 



a(t)a'(to) 



a 2 (t ] 



At 



H + z) 



a'(t ( 



a{t ] 



Ar. 



(2.12) 



where we used At = Ar for photons in the metric (2.2). Then the galaxies in the 
domain [r — \ Ar, r + | Ar] can be counted by 



5N = A7rp 9 f(r)Ar = An Pg f 



a(t( 



a' (t ) 



5z 
l + z : 



(2.13) 



where p g is the number density of galaxies in comoving volume dV = f ' 2 sin 9drd9d(f), 
which is a constant independent of t. By (2.13), we get the following rigorous equation 



[1 + z) 



dN_ 

dz 



a'(t )r 



a(t Q ) 



(2.14) 



By (2.9) and (2.11), we can solve r = r(t,z), then get the following function 

a'(t )r 



V = 



a(t ) 



(2.15) 



Generally speaking, ip = ip(t, z) is a complicated function. However for a small redshift, 
by the empirical Hubble's law, we have 



a'(t )r 




a' (t ) 


a(t ) 




a(t ) 2 



D r = H(t )D r =z, 



(2.16) 



where D r = a(to)r is the distance corresponding to radial coordinate r at time to- For 
flat space, the number of galaxies in a ball with radius r is given by 



Ar/ , 47T o 47T , 

N ( r ) = Y p 9 r = -JPgJ r - 



(2.17) 



Substituting (2.16) and (2.17) into (2.14), we get the distributive equation for flat space 

dN ■ 3 (2.18) 



Ndz z(l + z) 



The solution of (2.18) is given by 



N{z)^N 1 



l + z 



(2.19) 



By the precision of Hubble's law, we learn (2.19) is a good approximation for galaxy 
count at lest within the range z < 1. 

However, for the flat space K = 0, we can derive more accurate N(z) for large 
range z as follows. By (2.9), we always have solution to high accuracy 

a{t)=-t 2 , (for 5° < t < 60°), (2.20) 

t > 5° is because only the evolution after galaxies formed is concerned here, then we 
can omit the influence of ao- The condition t < 60° ~ l(rad) is derived by omitting 
the higher order term in the series (2.9). 
Substituting (2.20) into (2.11), we get 

\2 +1 



(t + r y-t 2 

t 2 



z = -^—f2 - = r (2 + r )> ( 2 - 21 ) 



or equivalently 



so we have 



ir^-^iwTTT (2 - 22) 

a(t ) t i + v ^TI 
Substituting (2.23) and (2.17) into (2.14), we get the distributive equation for flat space 
to high accuracy 



dN _ 3(1 + VTT^) 
~Wdz ~ 2z(l + z) 

The solution of (2.24) is given by 



(2.24) 



N{z) = N ° \TT7T7m) ■ (2 ' 25) 

Define the reference function by 

n{z) = In {^^yj = 2 HVTTz' - 1) - \ In(l + z) - n , (2.26) 

where z is the first redsheft point measured, from which we count galaxies, 

n = 21n( v / TT^- 1) - ^ln(l + z ) (2.27) 

is a constant. The function n(z) is the criterion which we are looking for, because for 
all z > zq, we definitely have 

5 

^(^)open > n(z) Rat = 2 ln(Vl + z - 1) - - ln(l + z) - n > n(z) doscd , (2.28) 
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The trend curves of n(z) 




Redshift z 



Figure 2: The Distributive Functions of Counting Galaxies 

due to the direct influence of f(r) on N(z). The trend curves of n(z) are displayed in 
Fig. (2). The practical counting data n(z) should also be a smooth function of z due to 
the law of large number. By this method, we will compare the trends of two smooth 
functions rather than two data, so there should be less ambiguous consequence caused 
by experiment and model errors. 

3 Discussion and conclusion 

From the above computation and analysis, we can get the following results. 

(CI). From Fig(l) we learn that the present mass energy Qk corresponding to the 
spatial type is generally small in value for a young universe, which is not sharp to give 
a definite answer for the spatial type. 

(C2). The space type can be determined by counting the galactic distribution with 
respect to the redshift z, which is a smooth function of z and can be sharply compared 
with the theoretical reference function (2.26). 

(C3). The cosmic equation of state p(a) is a useful concept for cosmology, because 
p is a superposable intensity with clear physical meanings, and it can be easily derived 
for concrete matter models via action principle, which seems to be more effective than 
w(z). 
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